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Abstract Current climate change scenarios predict 
that aquatic systems will experience increases in 
temperature, thermal stratification, water column 
stability and in some regions, greater precipitation. 
These factors have been associated with promoting 
cyanobacterial blooms. However, limited data exist on 
how cyanobacterial composition and toxin production 
will be affected. Using a shallow eutrophic lake, we 
investigated how precipitation intensity and extended 
droughts influenced: (1) physical and chemical condi- 
tions, (11) cyanobacterial community succession, and 
(111) toxin production by Microcystis. Moderate levels 
of nitrate related to intermittent high rainfall during 
the summer of 2013-2014, lead to the dominance of 
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Aphanizomenon gracile and Dolichospermum cras- 
sum (without heterocytes). Microcystis aeruginosa 
blooms occurred when ammonium concentrations and 
water temperature increased, and total nitrogen:total 
phosphorus ratios were low. In contrast, an extended 
drought (2014-2015 summer) resulted in prolonged 
stratification, increased dissolved reactive phospho- 
rus, and low dissolved inorganic nitrogen concentra- 
tions. All A. gracile and D. crassum filaments 
contained heterocytes, M. aeruginosa density 
remained low, and the picocyanobacteria Aphano- 
capsa was abundant. A positive relationship 
(P < 0.001) was identified between microcystin quo- 
tas and surface water temperature. These results 
highlight the complex successional interplay of 
cyanobacteria species and demonstrated the impor- 
tance of climate through its effect on nutrient concen- 
trations, water temperature, and stratification. 


Keywords Aphanizomenon - Dolichospermum - 
Microcystis - Surface water temperature - Nitrogen - 
Phosphorus 


Introduction 


The occurrence of cyanobacterial blooms is increasing 
in lakes across the globe (Garcia-Pichel et al., 2003; 
Paerl & Otten, 2013; Harke et al., 2016). Cyanobac- 
terial blooms can have negative environmental 
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impacts including reductions in water clarity and 
dissolved oxygen, which may suppress macrophyte 
growth and affect invertebrates and fish (Smith, 2003; 
Paerl, 2014). Many bloom-forming species also pro- 
duce natural toxins (cyanotoxins) which pose chronic 
and acute health risks to humans and animals (Chorus, 
2001; Codd et al., 2005). Blooms and their associated 
toxins also have economic and social consequences 
such as loss of tourism opportunities and increased 
costs associated with technical upgrades to water 
treatment plants and additional surveillance responsi- 
bilities (Kouzminov et al., 2007; Steffensen, 2008; 
Hamilton et al., 2014). 

Increases in anthropogenic eutrophication have his- 
torically been considered to be the main drivers of 
cyanobacterial proliferations, with the relative impor- 
tance of both nitrogen (N) and phosphorus (P) inputs 
debated for many decades (e.g. Smith, 1983; Lewis & 
Wurtsbaugh, 2008; Schindler et al., 2008; Conley et al., 
2009; Ahn et al., 2011; Smith et al., 2016). In recent 
years, the impacts of global climate change on the 
distribution and intensity of cyanobacterial blooms has 
attracted considerable scientific and media attention 
(Paerl & Huisman, 2008, 2009; Brookes & Carey, 2011; 
Carey et al., 2012; Kosten et al., 2012; Paerl & Paul, 
2012). Most climate change models predict that aquatic 
systems will experience increases in temperature, ther- 
mal stratification, and water column stability, all factors 
associated with promoting cyanobacterial blooms (Jöhnk 
et al., 2008; Wagner & Adrian, 2009; Trolle et al., 2010; 
Kosten et al., 2012). Less attention has been given to how 
other variables linked to climate change, such as 
variations in precipitation, extended droughts, or 
increased carbon dioxide levels, might affect cyanobac- 
terial biomass (O’Neil et al., 2012; Reichwaldt & 
Ghadouani, 2012). For example, prolonged and intense 
precipitation events can mobilise nutrients and sediments 
on land and increase nutrient enrichment of receiving 
waterbodies (Paerl et al., 2006). 

Understanding the effects of climate change on 
cyanobacterial blooms is complicated by successional 
patterns in cyanobacterial communities. For example, 
Fernandez et al. (2015) tracked the succession of 
Microcystis spp. and Dolichospermum (basionym An- 
abaena) circinale (Rabenhorst ex Bornet & Flahault) P. 
Wacklin, L. Hoffmann & J. Komárek during a bloom, 
and identified P, temperature and solar radiation as key 
factors in regulating shifts in their relative composition. 
Despite this and other examples (e.g. Rajaniemi- 
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Wacklin et al., 2008; Ni et al., 2012; Cai & Kong, 
2013), most research to date has focused on monospeci- 
fic blooms with the consequence that knowledge on 
cyanobacterial succession is limited (Fernández et al., 
2015). A further layer of complexity is that among 
cyanobacterial species known to produce cyanotoxins, 
both toxic and non-toxic genotypes can exist. Blooms 
are usually comprised both, and these can currently only 
be distinguished in the natural environment via molec- 
ular techniques (Kurmayer & Kutzenberger, 2003). 
Factors that promote one genotype to dominate over 
another are also largely unknown (Rinta-Kanto et al., 
2009; Van de Waal et al., 2011), and predicting how 
these might shift with climate change is challenging 
(Wood et al., 2015). Furthermore, the amount of toxin 
produced by toxic strains is variable, and factors such as 
temperature (Dziallas & Grossart, 2011; Kleinteich 
etal., 2012) and nutrients (Harke & Gobler, 2013; Horst 
et al., 2014) may be involved in regulating production 
and therefore also influenced by climate change. 

To date, most research investigating how variables 
linked with climate change might affect cyanobacte- 
rial biomass, species succession, genotype abundance, 
and toxin production has been based on models or 
culture-based studies, and there is uncertainty around 
the applicability of these data to real-world situations 
(Wood et al., 2015). Natural population surveys 
undertaken over consecutive years with contrasting 
and extreme climatic conditions may provide valuable 
insights and begin to address this uncertainty. Lake 
Rotorua (South Island, New Zealand; Fig. 1) is a 
shallow eutrophic lake which experiences annual 
cyanobacterial blooms between spring and autumn 
(Wood et al., 2011, 2012a). During the summer of 
2013-2014, the North Canterbury region of New 
Zealand (where Lake Rotorua is situated) experienced 
a series of rainfall events including intense precipita- 
tion from tropical cyclones Lusi (March 2014) and Ita 
(April 2014) as they transitioned into sub-tropical 
weather systems over New Zealand. In contrast, this 
region experienced a severe drought in the summer of 
2014-2015, with the second lowest total rainfall (over 
the summer period) measured since records began 
75 years before present. These contrasting weather 
patterns provided a unique opportunity to explore how 
differences in climate affect physical and chemical 
parameters in a shallow eutrophic lake, and the 
subsequent influence on cyanobacterial biomass, 
community composition, and toxin concentrations. 
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Fig. 1 Map of Lake 
Rotorua (South Island, New 
Zealand) showing the 
sampling site and thermistor 
chain sites 1 (January—April 
2014) and 2 (May 2014— 
May 2015). Inset map of 
New Zealand showing the 
location of the lake 


The aim of this study was to use a dataset collected 
from Lake Rotorua at weekly or fortnightly intervals 
between January 2014 and May 2015 to address the 
following three questions: (1) How are physical and 
chemical conditions within a shallow eutrophic lake 
affected by contrasting temperature and precipitation 
regimes Over consecutive summers? (11) Do physical and 
chemical changes affect cyanobacterial biomass, com- 
position, and succession, and which parameter(s) are the 
key drivers of any observed changes? (111) Did the relative 
abundance of toxic and non-toxic Microcystis genotypes, 
and toxin quota vary between the two summers and what 
were the possible drivers of these differences? Our final 
question focused specifically on Microcystis, as this is the 
only known toxin-producing cyanobacteria genus in 
Lake Rotorua (Wood et al., 2011, 2012a). 


Materials and methods 


Study site and sample collection 


Lake Rotorua (South Island, New Zealand 42°24’05S, 
173°34'57E) is a small (0.55 km”), Shallow (max. 
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Lake outflow 


173°35'E 


depth 3 m), eutrophic lake in the northeast of the 
South Island of New Zealand (Flint, 1975). The native 
podocarp-mixed hardwood forest that once dominated 
the catchment was removed in the early 1900s. The 
catchment now comprises mixed land use of low- 
intensity grazing, regenerating native scrub and bush- 
land, and exotic weeds and trees. Inflow is largely 
from surface run-off during rainfall events. During 
high rainfall, there is discharge from a small outflow at 
the southern end of the lake (Fig. 1). The changes in 
catchment land-use, resulting in greater nitrate and 
sediment run-off, coupled with a long residence time 
are likely reasons for the eutrophication of this lake. 
With its close proximity to the sea, Lake Rotorua 
experiences a relatively mild climate. The average 
temperature in the nearby town of Kaikoura is 12.7°C 
and precipitation averages 881 mm. The climate in 
this region is classified as Cfb by the Köppen-Geiger 
system. 

Temperature was measured at 5- or 15-min inter- 
vals using temperature loggers (HOBO®, Onset). 
Initially, this was undertaken from the end of a 
floating pontoon (January—April 2014; Fig. 1) with a 
chain of loggers at the lake surface (0 m), and 0.3, 
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0.45, 0.9, and 1.5 m below the water surface. There- 
after, the chain was positioned at a central point in the 
lake (Fig. 1) with loggers at the lake surface, and 0.8, 
1.0, 1.4, 1.7, 2.0, and 2.6 m below the surface. To 
enable the data loggers to be maintained at a constant 
distance from the surface regardless of fluctuations in 
lake level, the thermistor chain was attached to a buoy 
floating on the surface of the lake. Daily mean rainfall 
and maximum air temperature data at the Kaikoura 
weather station (situated approximately 8 km from the 
study site) were obtained from the Cliflow database of 
the National Institute of Water and Atmospheric 
Research, New Zealand. 

Between January 2014 and May 2015, surface 
water samples (400 ml) were collected weekly or 
fortnightly from the end of a floating pontoon (15 min 
length) on the eastern side of the lake (Fig. 1). 
Samples were stored chilled and in the dark, and were 
received at the laboratory within 24 h. Samples were 
subsequently well mixed and six sub-samples were 
taken: (1) 45 ml of pre-filtered lake water (20-um 
netting) was syringe-filtered (Whatman GF/C), and 
the filtrate was frozen (—20°C) for dissolved nutrients 
analysis; (2) 45 ml of lake water was frozen (—20°C) 
for total nutrient analysis; (3) 15 ml of lake water was 
preserved using Lugol’s iodine for microscopic anal- 
ysis; (4) 15—60 ml was filtered onto Whatman GF/C 
filters which were stored frozen (—20°C) for DNA 
extraction; (5) 5 ml of pre-filtered lake water (20-um 
netting) was syringe-filtered (Whatman GF/C), and 
the filtrate was frozen (—20°C) for extracellular 
microcystin analysis; and (6) 5 ml of lake water was 
frozen (—20°C) for total microcystin analysis. 


Laboratory analysis of weekly/fortnightly samples 


Total and filtered nutrients were analysed on a Lachat 
QuickChem® flow injection analyser (FIA+ 8000 
Series, Zellweger Analytics, Inc.) using the methods 
provided in APHA (2005). The nutrient species were 
analysed as ammonium (NH,-N), nitrite (NO.-N), 
nitrate (NO3-N), total N (TN), dissolved reactive P 
(DRP), and total P (TP). The limits of quantification 
were 0.002 mg 1~' for NH,-N, 0.001 mg 1~' NO,-N, 
and NO;-N, and 0.004 mg 1~' for TN, DRP, and TP. 

Cyanobacterial identification and enumeration was 
undertaken using an inverted microscope (1X70, 
Olympus). To assist in accurate enumeration, samples 
were mechanically ground (Tissue Grinder, Wheaton, 
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USA) for ca. 30 s to break up Microcystis colonies and 
filamentous species. Sub-samples (0.5-1 ml) were 
settled in Utermöhl chambers (1958), and cyanobac- 
teria cells and heterocytes were enumerated by 
scanning one-two transects at x400—600 magnifica- 
tion. Species identification was undertaken with 
reference to Baker and Fabbro (2002), McGregor & 
Fabbro (2001), Wood et al. (2005), and Komárek & 
Anagnostidis (1999). A second sub-sample (un- 
ground) from approximately every second sample 
was prepared as described above, and all non- 
cyanobacterial phytoplankton were identified and 
enumerated. Biovolumes were determined for each 
species using median cell dimensions (n = 30-50) 
measured at x 1,000 magnification under oil immer- 
sion (Olympus BX51) and volumetric equations of 
geometric shapes closest to each cell shape (Sun & 
Liu, 2003). 

Microcystis sp. strain CAWBG617, isolated from 
Lake Rotorua (Rogers, 2014), was grown in a glass 
flask (500 ml) in MLA medium (Bolch & Blackburn, 
1996) under a light regime of 90 umol s~' m~? with a 
12:12-h light/dark cycle, at a temperature of 
18 + 1°C. When the culture was in the exponential 
growth phase, a sub-sample (60 ml) was filtered 
(Whatman GF/C) and the filter stored frozen 
(—20°C) for DNA extraction. This was used to prepare 
a standard curve for the quantitative polymerase chain 
reaction (QPCR) assay used to assess toxic genotype 
concentrations. A second sub-sample (10 ml) was 
preserved using Lugol’s iodine and the cell concen- 
tration determined by microscopy (as described 
above). 

DNA was extracted from GF/C filters using a power 
soil DNA isolation kit (Mo Bio, USA) according to the 
protocol supplied by the manufacturer. All DNA 
samples from the Lake Rotorua samples were 
screened in duplicate for inhibition using an internal 
control assay. Each 12.5 ul reaction contained 6.25 ul 
Kapa probe fast QPCR kit master mix (x2), 1 ul of 
primers targeting the internal transcribed spacer 
region 2 of the rRNA gene operon of Oncorhynchus 
keta salmon sperm (0.4 uM, Sketa F2 and Sketa R3, 
Integrated DNA Technologies, USA Haugland et al., 
2005), 1 ul TaqMan probe synthesised with a FAM 
reporter dye at the 5’-end and a Black Hole Quencher 2 
at the 3’-end (0.2 uM, Sketa P2, Integrated DNA 
Technologies, USA; Haugland et al., 2005), 1 ul 
extracted salmon sperm DNA (15 ng, Sigma, USA), 
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and 1 ul of template DNA. The cycling profile was 
95°C for 3 min, followed by 50 cycles at 95°C for 3 s 
and 58°C for 10 s. When inhibition was observed 
samples were diluted (1/10) and re-analysed. 

QPCR, used to enumerate to the copy numbers of 
mcyE genes, was undertaken in triplicate for each 
sample in a 12.5 ul of reaction mix containing 6.25 ul 
Kapa probe fast QPCR kit master mix (x2), 1 pl of 
primers targeting a region within the mcyE open 
reading frame of the microcystin synthase gene 
(0.4 uM, mcyE-F2 and MicmcyE-R8; Vaitomaa 
et al., 2003), 0.2 ul of McyE probe (Rueckert & Cary, 
2009), and 1 ul of template DNA per sample. DNA 
from CAWBG617 was used to generate five-point 
standard curves ranging from 10.8 x 10° to 
10.8 x 10° cells ml~’. Each point of the standard 
curve was analysed in triplicate for each QPCR run 
conducted. The standard curve generated was linear 
(R? > 0.99) and PCR efficiency was >0.8. 

Sub-samples collected for extracellular microcystin 
were centrifuged (10,000 g, 5 min), and the super- 
natant was placed in a septum-capped vial. Sub- 
samples for total microcystin analysis were freeze- 
thawed and sonicated (30 min, 60 kHz) four times 
after adding formic acid (final concentration 0.1% 
v/v). The extract was clarified by centrifugation 
(10,000 x g, 5 min) and the supernatant was placed in 
a septum-capped vial. Microcystin samples were 
analysed directly or diluted (1/10-1/100 with 50% 
methanol containing 0.1% formic acid) by liquid 
chromatography—tandem mass spectrometry as 
described in Wood et al. (2012b) and Clearwater 
et al. (2014). Primary standards of the microcystin 
congeners -RR, -YR, and -LR were purchased from 
DHI Lab Products (Denmark). The amount of micro- 
cystin per toxic cell was calculated by summing the 
concentration of all congeners in the samples, sub- 
tracting the extracellular microcystin values from the 
total microcystin values, and dividing this value by the 
mcyE copies ml~' (as determined using QPCR). 


Data and statistical analysis 


Thermocline depth was quantified using the rLakeA- 
nalyzer package with the software R with a daily 
smoothing function applied (R Development Core 
Team, 2014). Differences in environmental variables 
between the summers of 2014 (for this study defined as 
January—May 2014) and 2015 (January—May 2015) 
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were evaluated. These were tested using a one-way 
ANOVA analyses with year as a fixed factor. 
Response variables were transformed if necessary, to 
fulfil the assumptions of normality and heterogeneity 
of variance. 

The relationship between cyanobacterial commu- 
nity composition (measured as biovolumes) and 
environmental variables was analysed using multi- 
variate multiple regression (McArdle & Anderson, 
2001), using the DistLM routine in the PRIMER 6 and 
PERMANOVA (Anderson et al., 2008). A marginal 
test was used where individual variables were fitted 
separately to test their relationship with the cyanobac- 
terial data (ignoring other variables), followed by a 
step-wise selection procedure, conditional on vari- 
ables already included in the model and using the 
Akaike information criterion (AIC) with a correction 
for finite sample sizes (AIC,) selection criteria. The 
conditional test identifies the subset of variables that 
best predicts the observed pattern in cyanobacterial 
community composition. Both the conditional and 
marginal tests were undertaken with 4,999 permuta- 
tions using Bray—Curtis similarities of the log + 1 
transformed cyanobacterial biovolume data. Drafts- 
man plots were used to check multi-skewness among 
environmental variables, and all were log + 1 trans- 
formed prior to the analysis and included: NO,-N, 
NO3-N, NHy-N, DRP, DIN:DRP, TN, TP, TN:TP, 
weekly averaged surface water temperature and total 
weekly rainfall. 

The relationship between Microcystis biovolume, 
percentage of toxic Microcystis cells, microcystin 
quota, and selected environmental variables was 
modelled using generalised linear models within the 
software R (Zuur et al., 2013; R Development Core 
Team, 2014). Microcystis biovolume and microcystin 
quotas were modelled using gamma errors and a log 
link, whereas the percentage of toxic Microcystis cells 
was modelled using beta error and a logit link 
(Pinheiro & Bates, 2000; Zuur et al., 2013). 

Initial data exploration was conducted to identify 
outliers, homogeneity of variances, normality, excess 
of zeros, collinearity, and the relationship between 
predictors and response variables following the pro- 
tocol described in Zuur et al. (2010). Covariates were 
transformed if necessary to avoid right-skewness and 
achieve unimodal distributions. Collinearity among 
predictor variables was initially checked by calculat- 
ing the variance inflation factor (VIF) for each 
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covariate and sequentially dropping the covariate with 
the highest VIF until all were <5 (Zuur et al., 2010). 
To avoid over-parametrization of the models only 
TN:TP, log(DIN:DRP) and surface water temperature 
were included. Additionally for the analysis of the 
percentage of toxic Microcystis cells and microcystin 
quota, Microcystis biovolume was included as a 
covariate in the initial models. Models were selected 
using a step-wise process based on generalised AIC 
and validated by inspecting the normalised quantile 
residuals. 


Results 


Surface water temperature in Lake Rotorua ranged 
from 5.7°C (24 July 2014) to 33°C (4 February 2014), 
and from 5.7°C (24 July 2014) to 25.9°C (22 February 
2014) at 2.6 m depth (Fig. 2). When the thermistor 
chain was positioned close to the shore (January—April 
2014) diurnal stratification was apparent, with the 
presence of a shallow surface thermocline (0.1—0.5 m; 
Fig. 2). From May 2014 to June 2015, when the 
thermistor chain was positioned mid-lake, a relatively 
persistent thermocline was present at ca. 1—-1.2 m from 
October to December 2014, and April-May 2015 
(Fig. 2), and at ca. 1.8-2 m in summer (January—May). 

From January and May 2014, total rainfall was 
540 mm compared to 82 mm during the same period 


Depth (m) 





15Jan14 22Feb14 O1Apr14 22 Jun 14 


Fig. 2 Temperature contour plot from Lake Rotorua (Kaik- 
oura, New Zealand) over the 17-month study period (January 
2014—May 2015). From January to April 2014, the thermistor 
chain was positioned off the end of the floating pontoon, and 
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in 2015. During this period, in 2014, there were five 
large rainfall events (>20 mm day’). The two largest 
were associated with tropical cyclones Lusi (15 
March, 79 mm) and Ita (17 April, 138 mm; Fig. 3). 

From January 2014 to May 2015, NH4-N concen- 
trations in the surface water were mostly low 
(<0.09 mg 1~'; Fig. 4A). Several exceptions occurred 
during the warmer water periods, on 26 February 2014 
(0.23 mg 17"), 30 April 2014 (0.11 mg 17'), 6 
November 2014 (0.4 mg 17'), 3 March 2015 (0.2 mg 
1~'), and 1 May 2015 (0.2 mg 17‘). During the winter 
period of July to August 2014, concentrations were 
generally above 0.2 mg 1"! (Fig. 4A). Nitrite and 
nitrate concentrations were low (<0.05 mg 1~') during 
both summer periods but increased markedly from 
March to September 2014 (0.1-0.7 mg 175. Surface 
water DRP concentrations were low (average 
0.005 mg 1~') from January 2014 to November 
2014, with the highest concentrations (0.02 mg I=) 
on 11 March 2014 and 6 November 2014 (Fig. 4B). 
Concentrations of DRP were higher from December 
2014 to May 2015 (average 0.02 mg 1~'), with three 
notable peaks on 29 December 2014 (0.05 mg 175,20 
January 2014 (0.04 mg 17'), and 26 March 2015 
(0.05 mg 17’). 

TP and TN were higher during the summer months 
(January—May) of both years (Fig. 4A, B). Over the 
period from January to May 2014, the DIN:DRP ratio 
was highly variable (average of 33), and in 2015, it 


Temperature (°C) 





15Oct14 30Jan15 09 Mar15 19 May 15 


thereafter, it was position at a central point in the lake (See 
Fig. 1 for locations). The white section on the plot is due to the 
shallower depth of the initial deployment. The black line shows 
the position of the thermocline 
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Fig. 3 Total daily rainfall at the Kaikoura Weather Station between January 2014 and May 2015. The red boxes show the two time 


periods used for comparative purposes in this study 


was consistently low (average of 6; Fig. 4C). The 
TN:TP ratio was <5 from January 2014 to April 2015 
and increased to >15 from May to July 2014. From 
November 2014 to the conclusion of the study, 
the TN:TP ratio was relatively constant between 10 
and 15 (Fig. 4C). 

The ANOVA analysis identified significant differ- 
ences in summer periods of 2014 and 2015 (assessing 
only the months of January—May for each year) for a 
number of physical and chemical parameters (Fig. 5). 
Higher surface water temperatures (P < 0.01), air 
temperature (P < 0.05), and lower rainfall (P < 0.05) 
were observed in 2015 compared to 2014 (Fig. 5). 
Comparing these two periods, there was a significant 
increase in the mean TN and DRP concentrations 
(P <0.01) and TN:TP ratio (P < 0.001), and a 
significant decrease in the DIN:DRP ratio (P < 0.01). 

Twenty-nine phytoplankton species were recorded 
in the lake; 3 bacillariophytes, 14 chlorophytes, 1 
crytophyta, | charophyta, 6 cyanobacteria, 2 dino- 
phytes, 1 euglenophyte, and 1 ochrophyta. Throughout 
the study period, the dominant phytoplankton phyla 
were varied; cyanobacteria were dominant at the start 
of the study until late February 2014; thereafter, 
chlorophytes (predominantly Dictyosphaerium sp. and 
Senedesmus spp.) or an euglenophyte (Trachelomonas 
sp.) were most abundant (Fig. 6). 

Of the six cyanobacterial species detected in Lake 
Rotorua, only four were detected in notable abun- 
dance: Aphanocapsa sp., Aphanizomenon gracile 
Lemmermann, Dolichospermum crassum (Lemmer- 
mann) P. Wacklin, L. Hoffmann & J. Komárek, and 


Microcystis aeruginosa (KĶützing) Kützing (Fig. 7A). 
Between 29 January and 22 February 2014, D. 
crassum dominated the waterbody, reaching a maxi- 
mum biovolume of 87 mm? 17' on 5 February 2014 
(Fig. 7A). A. gracile was present at moderate levels 
(<8 mm? 17!) and Aphanocapsa sp. at low levels 
(<0.1 mm? 1~') during this period (Fig. 7A). Micro- 
cystis was observed in samples from 5 February 2014 
and became dominant on 5 March 2014 (Fig. 7A). 

Over the winter—early spring period of June- 
October 2014, total cyanobacterial biovolumes were 
<1 mm? 17! (Fig. 7A). From 15 October 2014 to 9 
February 2015, the majority of the total cyanobacterial 
biovolume was A. gracile (Fig. 7A). From 25 Febru- 
ary 2015, D. crassum accounted for the greatest 
biovolume of the four species for the remainder of the 
study period, reaching a maximum of 23 mm? 17! on 
22 April 2015. This was in contrast to the same period 
of the previous summer when D. crassum was only 
present at low levels (<1 mm? 17'; Fig. 7A). M. 
aeruginosa Was present in low concentrations 
(<0.8 mm? 1~') from 23 October to 26 November 
2014. It was not present again until 29 December 2014 
and was then observed in all samples until the end of 
the study, but M. aeruginosa biovolumes never 
exceeded 3 mm? 17! (Fig. 7A). Between 9 February 
and 26 March 2015, Aphanocapsa was present in high 
concentrations, with biovolumes of ca. 1 mm? 17! and 
a maximum of 3.5 mm? 17! (Fig. 7A). 

Between January and May 2014, heterocytes 
occurred in low abundance in D. crassum filaments 
and none were observed in A. gracile filaments 
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Fig. 4 Nutrient concentrations for Lake Rotorua (Kaikoura, 
New Zealand) over the 17-month study period of January 2014— 
May 2015; A ammonium (NH,-N), nitrate + nitrite (NO2 + 
NO;-N), and total nitrogen (TN) concentrations, B dissolved 


(Fig. 7B). In contrast, heterocytes were common in D. 
crassum and A. gracile filaments between October 
2014 and May 2015 (Fig. 7B). Interestingly, hetero- 
cyte frequency in D. crassum was highest prior to the 
peak biomass of this species, whereas the highest 
heterocyte frequency in A. gracile occurred after its 
peak biomass (Fig. 7B). 

Microcystins were detected in 43% of the samples 
(Fig. 7C). Six microcystin variants were detected 
during the course of the study: -RR, desmethyl-LR, 
didesmethyl-LR, -FR, -WR, and -LA, although -FR, - 
WR and -LA were only detected on one occasion (26 
April 2014; data not shown). Total microcystin 
concentrations were generally higher in 2014 than 
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reactive phosphorus (DRP) and total phosphorus (TP) concen- 
trations, and C TN:TP and dissolved inorganic nitrogen 
(DIN):DRP 


2015, with the maximum concentration of 7.2 ug [= 
measured in the sample from 19 March 2014 
(Fig. 7C). The highest microcystin quotas were mea- 
sured between 29 January 2014 and 12 February 2014 
(1.14-2.14 pg cell”'; Fig. 7C) with all other quotas 
<0.54 pg cell~'. Marked temporal variability in 
microcystin quotas was observed, for example, there 
was a 280-fold decrease within a 2-week period 
(2.14 pg cell”! on 12 February 2014 to 0.01 pg cell™! 
26 February 2014; Fig. 7C). 

Microcystis genotype composition was highly 
variable in 2014 (between 1 and 68% were toxic; 
Fig. 7D). Although no microcystins were detected, 
toxic genotypes were measured using QPCR between 
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Fig. 5 Boxplot of environmental variables between January 
and May for 2014 and 2015. Solid black line shows median, box 
shows first and third quartiles, whiskers extend to the last data 
point within 1.5 times the inter-quartile range where those data 
exist. Black dots are outliers beyond this range. Note the log 
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Fig. 6 Phytoplankton biovolumes in water samples from Lake Rotorua (Kaikoura, New Zealand) over the 17-month study period 
(January 2014—May 2015). Other Bacillariophyta, Dinophyta, Ochrophyta, Charophyta, and Cryptophyta 


29 October and 26 November 2015 (1-62%; Fig. 7D). 
During the summer of 2015, the abundance of toxic 
genotypes was more stable (32-76%) declining in 
relative abundance from the end of March 2015 to the 
conclusion of the study (Fig. 7D). 


Multivariate regression analysis identified six 
parameters that individually had a significant relation- 
ship with the cyanobacterial community composition: 
NO3-N, TN, TP, TN:TP, DIN:DRP and surface water 
temperature (Table 1A). The sequential model 
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Fig. 7 A Cyanobacterial biovolumes, B heterocyte frequency 
(HF), C microcystin (MC) quotas of toxic genotypes and total 
microcystin concentrations, and D Microcystis concentrations 


showed that the four variables together explained 
53.6% of the total variation of cyanobacterial com- 
munity composition (Table 1). The univariate regres- 
sion analyses indicated a significant (P < 0.001) 
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2015 


and percentage toxic genotypes in water samples from Lake 
Rotorua (Kaikoura, New Zealand) over the 17-month study 
period (January 2014—May 2015) 


negative relationship between Microcystis biovolume 
and TN:TP (Fig. 8A). 

Univariate analysis identified a significant positive 
relationship (P < 0.001) between surface water 
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temperature and microcystin quota (Fig. 8B). There 
were no significant relationships between the environ- 
mental variables explored and the percentage of toxic 
Microcystis cells. 


Discussion 


Variability in physical and chemical conditions 
between summers 


An increase in the frequency and intensity of rainfall 
events is one of the predicted outcomes of climate 
change for some regions of the globe (IPCC, 2007). 
However, few studies have investigated the direct 
effect of rainfall on cyanobacterial biomass and 
composition (reviewed in Reichwaldt & Ghadouani, 
2012). High intensity rainfall events have been shown 
to cause either a decrease in biomass, or a complete 
collapse of the bloom due to de-stratification, increased 
turbidity, dilution, and flushing (Jacobsen & Simonsen, 
1993; Jones & Poplawski, 1998). These studies concur 
with our observations in Lake Rotorua, where the most 


Table 1 Results of non-parametric multiple regression of 
cyanobacterial community composition on individual environ- 
mental variables for: (A) each variable taken individually 


Variables F 


(A) Marginal test 


Nitrite 1.3 
Nitrate 10.9 
Dissolved reactive phosphorus (DRP) 1.7 
Ammonium 0.9 
Total nitrogen (TN) 21.6 
Total phosphorus (TP) 23.2 
TN:TP 5.0 
Dissolved inorganic nitrogen:DRP 4.9 
Surface water temperature 10.9 
Rain 0.6 
(B) Sequential test 
TP 23.2 
Temperature 13.6 
TN 5.0 
DRP 3.4 
Nitrite 22 


P values in bold are considered significant (<0.05) 
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extreme rainfall event (17 April 2014) led to cooling 
and complete mixing of the water column (Fig. 2) and 
coincided with a rapid decrease in cyanobacterial 
biomass (<3 mm? 1~!, which was low on subsequent 
sampling occasions for around 6 months, until spring; 
Fig. 7A). Although discharge in the lake outflow was 
not directly measured, the authors were at the lake on 
18 April 2014 and noted an increase in lake level by 
>] m from observations on the previous day and 
substantial discharge in the previously dry outflow. 
Both the lake water outflows had become highly turbid. 
In contrast, over the same period in 2015 cyanobac- 
terial biovolume remained high (>22 mm? 1~') until 
sampling ceased in June (Fig. 7A). 

Both prolonged rainfall and intense rainfall events 
may promote cyanobacterial blooms through increased 
particulate and soluble nutrient loads entering lakes 
(Toth & Padisak, 1986; Bormans et al., 2005). In Lake 
Rotorua, there was a marked spike in NO, + NO3-N 
concentrations in the first sampling (30 April 2014) 
which was conducted after the 17 April 2014 rainfall 
event (Figs. 3, 4A). Changes in TP were not observed, 
suggesting that most of the particulate material flushed 


(ignoring other variables) and (B) step-wise selection of 
variables, where variation explained by each variable is added 
to the model conditional on variables already in the model 


P % Var % Cum. 
0.28 2.4 
<0.001 17.3 
0.16 3.3 
0.43 1.7 
<0.001 29.3 
<0.001 30.8 
<0.01 8.7 
<0.01 8.6 
<0.001 17.4 
0.63 1.1 
<0.001 30.8 30.8 
<0.001 14.5 45.4 
<0.01 5.0 50.4 
<0.05 3.3 53.6 
0.08 2.1 55.8 


% Var percentage of variance in species data explained by that variable, %Cum. cumulative percentage of variance explained 
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Fig. 8 Relationships between A Microcystis biovolume and ratio of total nitrogen—total phosphorus (TN:TP), and B microcystin (MC) 
quota and surface water temperature. The shaded area shows the pointwise 95% confidence interval of the fitted values 


into the system had settled within this 2-week period. 
Interestingly, possible effects of run-off from less 
intense rainfall events earlier in the summer were not 
apparent in the lake nutrient analysis. For example, 
there was no increase in NO2 + NO3-N concentrations 
following a significant (79 mm) rainfall event on 15 
March 2014 (Figs. 3, 4A). We speculate that the high 
cyanobacterial and phytoplankton biomass in the lake 
at this time would have rapidly taken up incoming 
soluble nutrients (Takamura et al., 1987), resulting in 
no detectable change in NO, + NO3-N in samples 
from 19 March 2014 (Fig. 4A). In winter, when 
cyanobacterial and other phytoplankton biomass was 
low, rainfall events (e.g. 9 June 2014) usually resulted 
in a marked increase in water column NO, + NO;3-N 
concentrations (Fig. 4A). Different forms and concen- 
trations of inorganic N may favour different algal 
groups. Blomqvist et al. (1994) suggested eukaryotic 
algae are favoured by NO3-N. This aligns with the 
observation of high abundances of chlorophytes and 
euglenophyte during the winter of 2014 when NO, + - 
NO3-N concentrations were elevated. Despite the 
decline in NO2 + NO3-N concentrations from Septem- 
ber 2014 to the end of the study, the biovolume of 
eukaryotic algae remained higher than cyanobacteria. 
This may partly be due to the significant increase in the 
TN:TP ratio over this period. The optimal ratio for 
eukaryotic algae is believed to be between 16 and 23 
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molecules of N:1 molecule of P, with ratios of 10-16:1 
favouring bloom-forming cyanobacteria (Mur et al., 
1999). 

One of the most notable differences in the cyanobac- 
terial community between the summer periods of 2014 
and 2015 was the almost complete absence of hetero- 
cytes in A. gracile and D. crassum filaments in 2014, 
whereas the heterocyte frequency (i.e. relative to vege- 
tative cells) was relatively high in these species 
throughout 2015 bloom (Fig. 7B). The development of 
heterocytes in response to low concentrations of dis- 
solved inorganic nitrogen (DIN) has been observed in the 
field (Laamanen & Kuosa, 2005; Walve & Larsson, 
2007; Wood et al., 2010) and in laboratory cultures 
(Zapomelova et al., 2008). The absence of substantial 
rainfall events in summer 2014 likely reduced the 
delivery of DIN from catchment sources and created 
strong DIN limitation on biomass, reinforcing domi- 
nance of N-fixing cyanobacteria during the bloom. 

Predicted increases in dry periods and increased air 
temperatures may promote the proliferation of many 
cyanobacterial species through higher water temper- 
atures and greater water column stability (Johnk et al., 
2008; Wagner & Adrian, 2009; Kosten et al., 2012). 
As predicted, the dry summer of 2015 resulted in 
significantly warmer air and surface water tempera- 
tures (Fig. 5). The two different positions of the 
thermistor chain in Lake Rotorua in 2014 and 2015 
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means a direct comparison of water column temper- 
atures were not possible. There was a prolonged and 
relatively stable period of stratification (ca. 1.9m 
depth) from mid-January to mid-March in 2015 
(Fig. 2). While we did not measure dissolved oxygen, 
we suggest that this stable stratification would have 
allowed development of an anoxic hypolimnia pro- 
moting the release of nutrients (e.g. DRP and NH,-N) 
from the bottom sediment (Chorus & Mur, 1999). 
Surface water concentrations of DRP were signifi- 
cantly higher in 2015, with the spikes in surface 
concentrations largely corresponding to periodic de- 
Stratification events (Fig. 2), likely caused by wind- 
driven mixing and convective cooling. 


Cyanobacterial biomass, composition, 
and succession 


In contrast to predictions from climate models and 
previous culture-/field-based studies, the warmer water 
temperatures, prolonged stratification, and higher 
concentrations of DRP in 2015 did not result in an 
increase in cyanobacterial biomass in Lake Rotorua. 
The average cyanobacterial biomass between January 
and May 2015 (10.7 mm? 1~') was almost one-half of 
that of the same period in 2014 (20 mm? 1~'). Closer 
examination of the cyanobacterial data provides sev- 
eral possible explanations to this conundrum. First, an 
average cyanobacterial biovolume of over 10 mm? 17! 
is still extremely high, and there was a persistent dense 
cyanobacterial bloom during most of the summer in 
2015. Second, the 2014 data are skewed by one 
extremely high biovolume value (87 mm? 17! or 
1.2 x 10° cells ml~'; 2 February 2014). Removal of 
this data point reduced the average biovolume across 
the January—May period to 15 mm? 1~'; however, as 
noted below this event appears to have had a profound 
influence on cyanobacterial succession in 2014. 

The multiple regression sequential model identified 
temperature, TP, and TN and as the three most 
important variables related to cyanobacteria compo- 
sition over the study period. The temperature at which 
maximum replication rate occurs varies among 
cyanobacterial species (Reynolds, 1989, 2006). Addi- 
tionally, changes in temperature also have indirect 
effects, or are a response to other events such as 
stratification, or cooling due to large rainfall events. 
There is little doubt that temperature plays a direct and 
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indirect role in structuring the cyanobacterial com- 
munity in Lake Rotorua, and this is discussed in 
further detail below for specific species. The impor- 
tance of N versus P in explaining total cyanobacterial 
biovolume and that of specific cyanobacterial taxa has 
been demonstrated and discussed for many decades 
(e.g. Smith, 1983; Watson et al., 1997; Vrede et al., 
2009; Dolman et al., 2012; Smith et al., 2016). 
Although the multiple regression identified N and P as 
important, careful consideration of temporal variabil- 
ity, internal, and external fluxes and uptake rates is 
needed to understand their role in structuring the 
succession of cyanobacterial communities. 

Within the two summers, there was a distinct 
succession in the cyanobacterial species. In 2014, D. 
crassum, and to a lesser extent A. gracile, initially 
dominated the biomass (Fig. 7A). Although both 
species are known N-fixers, the near complete absence 
of heterocytes in 2014 indicated that no active fixation 
was occurring and that there was sufficient DIN to 
support growth. In the sampling immediately after the 
extremely high biomass event on 2 February 2014 
(which was predominantly D. crassum), there was a 
marked increase in NH,-N, which may have resulted 
from the lysis of at least some of this biomass (as 
previously noted in Wood et al., 2010). Microcystis is 
extremely adept at utilising NHy4-N (Blomqvist et al., 
1994; Dai et al., 2012), and this in concert with other 
promoters of M. aeruginosa dominance such as high 
water temperature (Chu et al., 2007; Davis et al., 
2009), was likely to have initiated the bloom of this 
species. After a peak in biomass in February, M. 
aeruginosa declined as DIN decreased. The rainfall 
event of 15 March 2014 appeared to provide new 
soluble nutrients, promoting M. aeruginosa growth; 
until cyclone Lusi (17 April 2014) curtailed the bloom, 
most likely through washout of the biomass. M. 
aeruginosa biomass was markedly lower in 2015, but 
peaks also coincided with increased concentrations of 
NH. -N in the water. However, in these instances, we 
suggest that the increased concentrations were likely 
the result of NH,-N release from the sediment, as they 
occurred with short-term mixing events following a 
period of sustained stratification. The likely causes for 
the lower Microcystis biomass observed in 2015, 
compared to 2014, are that the 2015 NHy,-N pulses 
were of lower magnitude and did not occur during a 
period of elevated water temperature, and the TN:TP 
ratio was consistently high (see “Discussion” section). 
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The univariate analysis indicated an inverse relation- 
ship between Microcystis biomass and a TN:TP ratio 
(Fig. 8A). Although early research suggested that this 
non-N-fixing species would be favoured by higher 
TN:TP ratios, there is now increasing evidence that 
corroborates our findings of high biomass associated 
with low TN: TP (e.g. Marinho & de Oliveira e Azevedo, 
2007; Lehman et al., 2008; Dolman et al., 2012). For 
example, Takamura et al. (1992) investigated a Micro- 
cystis bloom in Lake Kasumigaura and showed that 
TN:TP ratios were mostly <10 during the bloom, but 
exceeded 20 after the disappearance of the bloom, a 
pattern also observed in Lake Rotorua. Beversdorf et al. 
(2013) noted a similar pattern to that observed in Lake 
Rotorua during their analysis of a 15-year dataset from 
the North Temperate Lakes Long Term Ecological 
Research Program. They found increased M. aerugi- 
nosa biomass associated with low N:P ratios, while 
Aphanizomenon, which potentially fixes N, was more 
likely to occur with higher N:P ratios. Marinho & de 
Oliveira e Azevedo (2007) provide evidence to show 
that low N:P ratios can be the consequence of high- 
uptake capacity and greater relative uptake rates of N 
than P by Microcystis, causing a decrease in the ratio. A 
similar result was also shown by Xie et al. (2003) 
although in their experimental setup both N and P were 
supplied in high concentrations. 

One of the most notable features of the cyanobacte- 
rial community in spring 2014 (October-December) 
was that it was comprised almost entirely of A. gracile. 
We did not sample over the same period in 2013, 
preventing an inter-annual comparison. During this 
period in 2014, the DIN:DRP ratio was extremely low, 
which is known to favour N-fixing species, and this may 
be one of the reasons for the dominance of A. gracile 
during this period (e.g. Vrede et al., 2009). Among 
N-fixing species, their capacity to fix N varies. For 
example, Dolman et al. (2012) noted differences in the 
N-fixation rate per heterocyte for a variety of nosto- 
calean species. A higher N-fixation capacity in concert 
with a preference for cooler temperatures (Reynolds, 
1989) could explain the dominance of A. gracile during 
this period. Allelopathy may be a further possible reason 
for its dominance as Borges et al. (2016) observed 
benthic Microcystis recruitment in Lake Rotorua was 
significantly inhibited when A. gracile was present in 
the sediment. Extracts from other Aphanizomenon 
species are known to inhibit cellular processes in other 
algae (e.g. Suikkanen et al., 2004). 
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Aphanocapsa sp. was present only in very low 
numbers in 2014, yet was abundant between January 
and April 2015. Aphanocapsa sp. is the smallest 
(diameter = 0.9 um) of the cyanobacteria species 
observed in Lake Rotorua. The large surface area— 
volume ratio of picocyanobacteria has been proposed 
as a reason for their dominance when nutrient 
concentrations are relatively low (Legendre & Ras- 
soulzadegan, 1995; Raven, 1998). Thus pico- 
cyanobacteria dominance in Lake Rotorua is 
consistent with low DIN concentrations between 
January and April 2015. 

Havens et al. (2016) examined the effect of extreme 
weather events and cyclical changes in water level 
using long-term (15-25 years) datasets from three 
shallow lakes: Lakes Okeechobee and George (Flor- 
ida, USA) and Lake Taihu (China). Tropical cyclones 
resulted in increases in run-off and sediment resus- 
pension, increases in dissolved nutrients, causing 
intense M. aeruginosa booms, findings that concur 
with many of our observations in Lake Rotorua. In 
contrast in Lake George, periods of high rainfall 
resulted in greater dissolved colour, reduced irradi- 
ance, and increased water turnover rates which 
suppress blooms, whereas in dry periods there were 
dense cyanobacteria blooms. Evaluating the degree to 
which extreme weather events affect shallow lakes is 
likely dependant on a myriad of interacting factors 
including: relative area to mean depths, catchment 
characteristics, sediment properties, and aquatic veg- 
etation. Additionally our data clearly shows that 
different cyanobacterial species have diverse 
responses to varying temperature and stratification, 
as well as different forms, concentrations and ratios of 
N:P. Our results point to the need to differentiate 
cyanobacterial biomass into species that better reflect 
their varied physiological attributes, as well to better 
resolve fluxes of nutrients that affect the competitive 
attributes of these species. 


Variability and drivers of toxic and non-toxic 
Microcystis genotypes and microcystin quotas 


Previous field-based studies have reported positive 
correlations between toxic Microcystis genotype 
abundance and TP (Rinta-Kanto et al., 2009; Yu 
et al., 2014), NO3-N (Yoshida et al., 2007) and water 
temperature (Conradie & Barnard, 2012). In the 
present study, no relationships were identified 


Hydrobiologia (2017) 785:71-89 


between the abundance of toxic genotypes and M. 
aeruginosa biomass or the environmental variables 
included in our model. To avoid over-parameterising, 
the model only a limited set of variables were 
incorporated, including those previously identified as 
important (or integrators of them such as TN:TP and 
DIN:DRP). Based on the data in the present study and 
findings from previous research that indicate a range 
of environment variables (e.g. temperature, nutrients) 
influence the relative abundance of toxic and non- 
toxic genotypes, we suggest that predicting these shifts 
in cyanobacterial blooms will be extremely challeng- 
ing. It is likely that variables that drive these changes 
vary temporally and spatially within and between 
lakes. 

During this study, marked variability in microcystin 
quota (>200-fold change) was observed over rela- 
tively short time frames (ca. 2 weeks). Some of this 
variation may be attributed to a community shift 
within the microcystin-producing strains as differ- 
ences can occur among toxic genotypes, with each 
producing different microcystin quotas and variants 
under the same conditions (Saker et al., 2005; Kriiger 
et al., 2010). As the QPCR assay used in this study 
only discriminates between toxic and non-toxic geno- 
types, no information on strain-level diversity was 
collected during the study. 

The univariate modelling used in the present study 
identified a significant positive relationship between 
microcystin quotas and surface water temperature. 
This corroborates earlier reports suggesting that 
microcystin concentrations are maximal around 
20—25°C (Van der Westhuizen & Eloff, 1985; Watan- 
abe & Oishi, 1985; van der Westhuizen et al., 1986; 
Amé & Wunderlin, 2005) and highlights the possibil- 
ity that predicted global increases in surface water 
temperature may lead to higher microcystin quotas. It 
is also possible that the higher surface water temper- 
atures could select for toxin-producing strains which 
produce higher microcystin quotas, but as noted above 
we were not able to ascertain this using the QPCR 
assay. 

The discrepancy between times of maximal Micro- 
cystis biovolumes, toxic genotype abundance, and 
microcystin quotas highlights the challenges for water 
managers trying to estimate periods of highest risk 
when relying solely on the cell concentration (or 
biovolume) data recommended by numerous guideli- 
nes (WHO, 2003; NHMRC, 2008; Wood et al., 2009). 
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For example, the second highest total microcystin 
concentrations detected in a water sample (3.6 pg 17’, 
30 April 2014) occurred when Microcystis biovolumes 
were moderate (2.8 mm? 1~'), the abundance of toxic 
genotypes was low (16%), and the microcystin quota 
was at intermediate levels (0.2 pg cell~'). To provide 
greater assurance in estimating risk from cyanotoxin 
exposure, parallel cyanobacterial identification, enu- 
meration, and toxin analysis are recommended. These 
data are essential to enable a progression from current 
models of Microcystis biomass or specific physiolog- 
ical attributes of cyanobacteria (see Oliver et al., 
2012), to in situ genotype composition and micro- 
cystin production, including how these variables 
might respond under predicted climatic change 
scenarios. 


Conclusions 


This study demonstrates the interactive effects that 
varying precipitation and temperature can have on 
physical, chemical, and biological parameters in 
shallow eutrophic lakes. The key differences observed 
in Lake Rotorua were higher NO3-N inputs during the 
wet summer, prolonged stratification, and increased 
DRP and low DIN concentrations in the drier warmer 
summer. Corresponding differences were observed in 
the cyanobacterial communities that reflected their 
varying physiological adaptations (e.g. N-fixation, 
nutrient uptake capacity, and temperature optima). No 
relationships between environmental parameters and 
shifts in toxic/non-toxic Microcystis genotypes were 
identified, although there was a positive relationship 
between temperature and microcystin quota. We 
suggest that while generalised prediction on how 
Microcystis biomass might alter under foreseen 
climate change scenarios are feasible, determining 
how genotype composition and microcystin produc- 
tion will vary is extremely challenging. 
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